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1. Introduction 

In stochastic mass transport processes, one is generally concerned with the dynamics 
of abstract, usually indistinguishable particles. Well known models include the totally 
asymmetric exclusion process (TASEP) [H |2] where the unidirectional free movement 
of particles is hindered by other particles, and the more generic zero-range process 
(ZRP) [3] that only models zero-ranged interactions of particles. Depending on the 
situation and given type of dynamics these particles may represent actual molecules on 
the microscopic scale, mesoscopic objects like intracellular motor proteins |1] or even 
macroscopic bodies such as people or cars in traffic systems |S]. Interesting effects such as 
phase transitions in the one-dimensional system and real-space condensation of particles 
can already be seen in the ZRP, despite the absence of ranged interactions ilElEli. 
In this paper we will consider a generalized transport process including short-ranged 
interactions which allows for a pair-factorized steady state (PFSS) |T0]. For a specihc 
type of interactions proposed by Waclaw et ah nn, this model features a spatially 
extended condensate that can assume several distinguished shapes. That is, by 
parametrization the condensate can be tuned to a single-site peak, rectangular or smooth 
parabolic envelope shape. Here we will consider and validate predictions on the scaling 
of the condensate width as well as its shape made in Refs. mm- 

The rest of the paper is organized as follows. In the next section we will 
introduce the considered stochastic transport model and summarize the derivation of 
the predictions. The methods needed to simulate the model and measure its properties 
will be discussed in the third section followed by a discussion of our results in the fourth 
section. The paper closes with our conclusions in the hfth section. 

2. Model 

Consider M indistinguishable particles on a one-dimensional, periodic lattice with N 
sites. Any site i can be occupied by any number of particles rui > 0. We only consider 
closed systems, that is, the number of particles is conserved, M = with an 

overall density p = M/N. The dynamics is dehned as a Markovian stochastic process: 
At any time step, a particle may leave from a randomly selected site and move to 
either of its neighbours. The specihc dynamics is largely controlled by the hopping rate 
u{mi\mi-i,mi+i), which determines whether a particle actually performs a hop. The 
fact, that it depends only on the occupation numbers of the selected site and those of 
its direct neighbours, rehects the short interaction range of particles and sites. 

The model can easily be tuned from symmetric to totally asymmetric dynamics by 
introducing a probability r, that is used to decide whether a particle hops to the right or 
left neighbour. Even though this allows tuning the process from equilibrium (r = 1/2) 
to non-equilibrium (r ^ 1/2) dynamics, the stationary state remains unaffected. 
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Such a model was proposed by Evans et al. HDl as an extension of the well known ZRP 
and therefore inherits several of its properties: The model features a stationary state 
that may contain a particle condensate. The steady state probabilities P{m) factorize 
over symmetric non-negative weight functions g{m,n) of pairs of sites 


P{m) = P(mi, . . . , TTlAr) 


1 

Z 


N 

Wg{mi,mi+i) 

i=l 


( 1 ) 


for = M kept constant, with the partition function Z = X]{m}-^(^) 

conhguration space {m}. In contrast, for the ZRP, the state probabilities factorize over 
single-site weight functions. 

If the weight function g{mi,mi 1) falls off faster than any power law, it has been 
shown that there exists a critical density p^, above which translational symmetry is 
broken and a particle condensate emerges [iniEiin]. That is, in the steady state the 
bulk system can only hold a limited number of particles given by and any further 
added particle increases only the condensate mass Me = M — pcN = (p — Pc)N■ 

The hopping rate that leads to the steady state Q is easily determined with the 
weight function g{m,n) as 




giruj - 1, mi_i) gimj - 1, m^+i) 

g{mi,mi-i) g{mi,nii+i) 


( 2 ) 


by resolving the balance condition in the steady state. In the case of symmetric hopping, 
detailed balance is fulhlled. 

An interesting way to construct such weights is to simplify and separate 


g{m,n) = \/p{m)p{n)K {\m — n\), 


(3) 


i.e., to factorize g{m,n) into a zero-range part p{m) and a short-range interaction part 
K{\m — n\) that depends on the difference of nearest-neighbour occupation numbers. 
The term zero-range interaction refers to the fact that it allows a particle on a given site 
to interact with particles on the given site only. In contrast, short-range interactions 
act between particles of different sites. For example, the zero-range process [3] with 
the hopping rate u{m) = 1 + h/m is, easily implemented using the weight function 
g{m,n) = = pim), with K{x) = 1 as there is no short-range interaction involved. 

In the original model by Evans et al. |T0] governed by the weights 


K{x) = exp (-J|x|), p{m) = exp (f/o<5m,o), 


(4) 


the condensate has a characteristic smooth parabolic shape. Under the conditions that 
K{x) falls off faster than any power law and p{m) = const for some m > the 

exact envelope shape can be calculated as well as its fluctuations and the condensate 
width scaling behaviour W oc y/M^ [12] • The factorization (§ of g{m,n) also makes 
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it easier to understand the physical picture of condensation involved. The zero-range 
potential p{m) gives a penalty to increasing occupation numbers, much like a potential 
energy term, while the short-range interaction K{x) tends to reduce the difference in 
occupation numbers of neighbouring sites, acting like a surface energy. It is then the 
respective relative importance of these terms that governs the properties of the emerging 
condensate. 

Waclaw et al. m then suggested a family of weights 

K{x) ~ e-“l^l^ p{m) ~ (5) 

that allows one to deliberately violate the asymptotic condition on p{m) for large 
m by choice of parameters. This produces regimes featuring qualitatively different 
condensation properties, and most notably may be used to tune the condensates envelope 
shape as well as the scaling behaviour of its width with system size [HI [12]. In Ref. [12] 
these regimes and properties are derived in the large-volume limit of the model. It is 
the subject of this work to check these predictions by means of numerical simulations 
at hnite volumes and scaling analyses. In the next subsection we will therefore shortly 
summarize the derivation of the interesting new properties. 

2.2. Theoretical predictions 

The proposed model is fully dehned by Eqs. 0.0 , and (|^. The parameters / 3,7 give 
the respective growth behaviour of the short-range and zero-range interaction strengths 
and are expected to have the most influence on condensation properties. The factors a, b 
allow to tune the relative strengths of the interactions but do not change the qualitative 
behaviour of the system as we will see later. In this paper we will therefore use a = b = 1. 

First, the condition 7 < 1 is derived to observe condensation at all. Next, the 
steady state weight Eq. Q can be thought to factor into one part accounting for the 
statistical weight of the condensate and one for that of the rest of the system. Then, 
a comparison of the weight of a single-site condensate Pi and that of a condensate 
occupying two or more sites shows that for /9 > 7 the weight grows faster than 
Pi, so that an extended condensate can be expected in the large system-size limit. 

In the case of the extended condensate a similar approach is used to predict some of 
its properties. Under the assumption, that the contribution of fluctuation of occupation 
numbers is small compared to the contribution from the condensate itself, the weight of 
a condensate extended over W sites is given as [H] 

InP(lU) lUc-l- ^InP ((mj+i — mi)) P ^lnp((mi)), ( 6 ) 

i i 

where for simplicity c = s — In Amax is treated as a positive constant to approximate 
the partition function, with an entropic term s that corresponds to fluctuations of the 
occupation numbers and the largest eigenvalue A^ax of the associated transfer matrix. 
The expectation values (ruj) correspond to values Hh{2i/W — 1) where h{f) with 
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f = 2i/W — 1 is the rescaled envelope shape of a condensate with height H = Mc/W 
and width W. Now, for the expectation of occnpation number differences (mj+i — mi) 
a distinction into two cases is proposed to allow simplifications. Either, the envelope 
shape h{f) is smooth or it is of rectangular shape. In the first case the expectation of 
the difference can be written using the derivative of the shape h'{f) and the expectation 
reads {2Mc/W^)h'{f). In the second suggested case, all differences are zero except at 
the condensate boundaries and the sum reduces to a single term. Then in both cases, 
the sums can be written as integrals using the rescaled shapes resulting in 


hlFsmooth(W") ~ VE 


df\nK 


\nPrect{W) ^Wc + 2\yiK 


deinp(Tr/i(0) , 

+W iyf\np{Hh{0)- 


(7) 

( 8 ) 


The sums decouple into a prefactor and a constant integral term, which does not depend 
on the condensate width. The prefactors however are used to determine the growth of 
the different condensate weights depending on their respective widths and masses 


InFsmoothW-VE 


/ 2 Mc 


c — a 


VhE2 


/3 rl 




(9) 


2 \h(l)\l>-bw(^\'l^dmor- ( 10 ) 

Finally, the condensate weights are estimatec||] by finding their maximal value with 
respect to the condensate width (with constant mass) which also gives the scaling laws 
for the condensate width 

In PsmoothiW) ~ lY _ ^ ^ 1 ^ 


In F,ect(H^) ~ W ~ for /3 > 0 


( 12 ) 


in the smooth and rectangular case, respectively. For 0 < /3 < 1/2, no solution exists for 
smooth condensates. Finally these weights are compared to find InFrect(hE) > InFsmooth 
for 1/2 < /d < 1, showing that the rectangular shape is the stable one in this regime. In 
Fig. 0 we plot the (Mc-independent) relative deviation 


Fdev (/3, y ) 


log(- log Fsmooth) - log(- log Prect) 
log(-logFrect) 


(13) 


to give an impression of the regime boundary. Figure shows these regimes and the 
respective condensate scaling laws. 


^peak 


0 , 


/3-7 

^smooth n ’ 

2/3-7 


Q^rect 


/3-7 

/3 -7 + 1 ’ 


(14) 


f During the estimation of the maximal weight for the given condensate type additional constraints 
are found where no solution exists, narrowing down the number of cases to look at. 
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Figure 1. Comparing the scaling behaviours of the most likely rectangular and smooth 
condensate shapes using the relative deviation (HI, over a region in /J-y-parameter 
space. Larger deviations i?dev (cut off at 0.06) mean that (in the thermodynamic 
limit) the probabilities are different, and lower ones that they are close. 
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Figure 2. The width-scaling exponent a in IF oc M“ for different parametrizations /3 
and 7 , cf. Eq. (14). The three regions of different scaling behaviour mark the distinctive 
phases. For 7 > 1 there is no condensation. This figure is a reproduction of Fig. 5 
in [T2] . 


3. Methods 

We will employ the empirical method of taking samples while evolving a Markov process. 
Once we have obtained a set of representative samples for one particular parametrization 
(/3, 7 , Me), we can determine the condensate width. With the widths for different 
condensate volumes Me, we then can estimate the condensate-width scaling exponent 

a. 
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3.1. Generating system samples 

The straightforward method for generating samples is to simnlate the process as dehned. 
Its dynamics with typically small transition probabilities leads to a large waiting time 
nntil the gathered samples are representative. The key ingredient for improving the 
simnlation is that one knows the steady state probabilities Q which are invariant 
nnder the employed dynamics, asymmetric hopping (r ^ 1/2) in non-equilibrinm or 
symmetric hopping (r = 1/2) in eqnilibrium. In the latter case, we can verify stationary 
eqnivalence to systems with other dynamics that yield the same distribntion, bnt mix 
mnch more qnickly. The Metropolis dynamics is known to mix well and can be readily 
adjnsted for the desired distribntion, which is why we chose to use it in our simulations. 
When the update scheme is altered so that particles may hop to any site instead of just 
the nearest neighbours, the simulation algorithm assumes the form: 

(i) Select randomly a site from which a move of one particle is attempted. 

(ii) If this site is empty, the system remains in the same state for another time step of 
1/N sweeps. If it is occupied, choose any of the other sites as a destination site. 

(iii) Accept the move with Metropolis probability Pace := min (Pfinai/Pinitiai, 1) or remain 
at the current state for another time step otherwise, where Pmitiai and Pfinai are the 
probabilities of the microstates before and after the proposed move, respectively. 

For our implementation of the simulation we use the reference implementation of a 
variant of the Mersenne Twister pseudo-random number generator [I3l[ll]. 

This approach works well for the smooth condensate shape region {/3 > 1), 
where the probability landscape is marginally rugged. For systems in the rectangular 
condensate regime this is not the case anymore because the probability landscape 
is heavily rugged as shown in Fig. The transitional states between rectangular 
condensates of different widths are suppressed by several orders of magnitude and 
simulations with single-particle (or local) updates mix slowly. 

In order to address this issue, we designed an additional set of multi-particle (non¬ 
local) updates that allow jumps between high probability regions, which correspond to 
the peaks in Fig. in a single step, thus bypassing the suppressed transitional states. 
This is achieved by defining a condensate expansion and reduction transformation, where 
the condensate width is increased or decreased by about one unit. 

With the above stated local update algorithm it is easy to verify that the resulting 
stationary distribution is given by Eqs. (§, (§, (g. By changing the dynamics we 
must be careful to keep the steady state unchanged. The local update performs well 
in exploring the landscape locally, which is why we will keep it and independently 
add non-local updates. The probability currents added by the new updates must 
cancel each other out at every state. This can be achieved by constructing a dynamics 
fulhlling detailed balance, where an incoming probability current is being canceled by 
its corresponding outgoing one. The one we used before and which we will use again is 
based on Metropolis transition rates. The main problem here is to treat equally the two 
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Figure 3. The shortest path in probability landscape starting from a condensate 
occupying a single site with 5000 particles and going to rectangular condensates of 
increasing widths. The peaks correspond to the high statistical weights of rectangular 
condensates and the valleys to those of the transitional states on a an ideal trajectory. 
The (unnormalized) probabilities are based on the parametrization with j5 = 0.6, 
7 = 0.4. 

non-local update moves and how often to use them in order to minimize the relaxation 
time. We chose to solve this by restricting the number of new transitions that are 
permitted at any state to (at most) two; an expansion and a corresponding reduction 
transition. The new simulation algorithm then consists of a preceding step, where a 
decision (according to some predefined probability) is made on whether to attempt a 
local or a non-local step[§] Choosing a local update is followed by the algorithm above, 
while a non-local update is governed by the following algorithm; 

(i) Decide whether to attempt an expansion or reduction transformation with equal 
probability. 

(ii) If the chosen transformation is not valid, the system remains at the current state 
for another time step (the system is constrained in a similar way to the situation 
of attempting to move a particle from an unoccupied site), otherwise continue. 

(iii) Accept the transformation with Metropolis probability Pace = min (PfinaiZ-Pinitiai, 1) 
or remain at the current state for another time step otherwise. 

From the algorithm above, the necessary conditions are that there should be only 
one way (if any) to expand or reduce the condensate and the transformation must be 
invertible such that reducing a previously expanded state always results in the original 
state. 

The conditions we decided on is that expansion transformations are eligible only if 
the original state has two unoccupied sites to the right of the condensate. One removes 
particles from the top of the condensate by sweeping from left to right and puts them 
on the empty site right to the condensate, until the two last condensate sites have the 
same occupation numbers (in some cases this is not possible, then the transformation 

§ For small probabilities it is better to use the waiting time, until a non-local update will be performed, 
which is geometrically distributed. 
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Figure 4. The transformation loop for an example condensate state including all 
intermediate steps. Any condensate state is suited for the respective transformation if 
it fits either starting point in this loop and adheres to the highlighted characteristics. 


is invalid). Likewise, reduction transformations are only valid if the original state has 
the same occupation numbers at the two rightmost condensate sites. One removes the 
particles at the rightmost condensate site and distributes them by sweeping from left to 
right onto the rest of the condensate. Figurej^gives an example of both transformations. 

For M = 1000 particles with fd = 0.7 and 7 = 0.4 the autocorrelation time of the 
condensate width of purely local updates is about 10 times larger than that of non¬ 
local updates. Increasing the particle volume by 50% (100%) increases this factor to 50 
(250). Whereas the new method enables systems with more than 10® particles to mix 
sufficiently in less than an hour, purely local updates are constrained to less than 10 ^ 
particles. 

In cases where fd and 7 are both close to 0, the critical density is so large, that the 
condensate is not surrounded by unoccupied sites anymore. Our non-local updates are 
then rendered ineffective and there is no easy way to adjust them to this situation. 

3.2. Determining the condensate width 

In the analytic approximation of the scaling behaviour ca. the width of the condensate 
is defined as the distance between the two outmost sites of the condensate, which we 
will call the condensate extension W (see Fig. |^a)). Whereas this is a good choice 
in the thermodynamic limit, for finite systems the condensate edges are subject to 
interactions with the rest of the system, which results in the edges being smeared 
outward the condensate. The condensate extension is an observable that is strongly 
influenced by these finite-size effects, and therefore overestimates the condensate width. 
This is especially the case with narrow but high condensates that have steep edges. 

We decided to measure the condensate width with another observable that is 
less prone to these smearing-out effects. This observable is based on the assumption 
that given a parametrization fd and 7 , the system exhibits a common characteristic 
shape that is merely rescaled in width and height in dependence on the volume of the 
condensate. Reversing this assumption, it must be possible to rescale condensate shapes 
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Figure 5. Comparison of mean condensate shapes at different volumes but the same 
parametrization /3 = 1.2 and 7 = 0.3. The shapes were rescaled in terms of extension 
by the condensate width and in terms of occupation numbers by normalization to 
a unit volume. With increasing condensate volumes the shapes converge to a single 
shape which indicates a common characteristic shape for this parametrization. The 
inset shows the original mean condensate shapes. 


(a) 



(b) 



Figure 6. Condensate width determination method by (a) taking the direct 
condensate base extension or (b) taking the distance of the left and right condensate 
wings respective centers of mass Wcm- 


of different volumes to such a common shape, which was done for one example system 
in Fig. 1^ One can see there, that as condensate volumes increase, the shapes converge 
to a characteristic shape, with small volume condensates deviating from it due to the 
mentioned hnite-size effects. 

First the condensate is separated in two wings left and right of its center of mass 
(CM) position icM = where ia and 4 are the begin and end positions 

of the condensate so that > 0 for i Then we calculate the width as the 

distance between the individual centers of mass of these two wings 


WcM 


Mr 


'I'b 

E 


irui 


i<iCM 




(15) 


Figure [^summarizes how this is done in comparison to the direct estimation method. 


3.3. Determining the condensate width scaling behaviour 

Once we have obtained estimates of the wing center-ov-mass width observable hFcM 
for system parametrizations with different condensate volumes Me, the prediction [ 12 ] 
states that for large volumes the scaling behaviour adheres to IFcm ot M" with q((/3, 7 ). 
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Figure 7. Log-log plot of the width scaling behaviour in dependence on the condensate 
volume Me for systems with /? = 0.9 and various 7 in the range [0.3,0.9]. The wing 
center-of-mass widths (15) have been normalized by dividing them by the smallest 
width. 



system size N 


Figure 8. The dependence of the condensate width with a fixed condensate volume 
Me on the system size for (3 = 1.3 and 7 = 0.5. The critical density is about pe — 0.14. 


Taking the logarithm this relation reduces to a linear relationship with a determining the 
slope. Figure [^conhrms that the data asymptotically approaches this linear relationship 
with increasing condensate volumes. Determining the scaling exponent is a linear 
regression problem, where the measured condensate width is uncertain but the measured 
condensate volume is considered certain. A minor problem is the overestimation of the 
condensate width for small condensates with strong smearing (hnite-size) effects, which 
we simply avoid by simulating large enough condensates. 

One potentially problematic aspect left untouched so far, is how the condensate 
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width and with it the scaling behavionr is affected by the system size, especially 
considering effective interactions between the condensate edges through the background. 
The measurements in Fig.j^show that for systems where the condensate volume stays the 
same, but the system size changes, the measured condensate width remains unchanged 
as long as p > Pc- This leads to the conclusion that for densities much larger than 
the critical density, the competition between the gaseous background and condensate 
phase is dominated by the latter. This system size independence is favorable for us, 
since it eliminates one degree of freedom in the discussion (leaving only the system 
parameterization /3, 7 and condensate volume M^) and allows us to choose systems 
sizes that are sufficient for the condensates, but small enough to maintain simulation 
performance. We used a fixed size about ten times the width of the largest condensate to 
minimize interaction of condensate edges around the boundary and eliminate percolating 
condensates. 


4. Results 

We are first going to look at the condensate phases, where especially the boundary 
between the smooth and rectangular condensate phase is of interest. Then we will 
present the estimated condensate width scaling behaviour. 


4 . 1 . Condensate phases 


In order to determine the phases we will use the previously introduced concept of 
characteristic shapes. Figure [^shows estimations of these shapes for a range of points in 
parameter space at a system volume of about 10^ masses. The data is obtained in runs 
with on the order of 10® Monte Carlo sweeps (consisting of N single-particle moves) 
taking 10^ samples of the system. For /3 < 1 we see characteristic shapes with constant 
elevation in the middle section and steep edges at the boundaries, which correspond 
to the predicted rectangular condensate shape. For small (3 and 7 we see that the 
estimation of the characteristic shape fails. This is due to the increased critical density 
and the associated failure of the non-local update. The detection method used here 
is unable to distinguish between extended rectangular condensates and peak-shaped 
single-site condensates. 

For /3 > 1 we see characteristic shapes with decreasingly steep edges that converge 
to a dome. These bell-shaped condensates correspond to the smooth case that was 
predicted for this region. For small 7 and /3 ~ 1 the condensation process is close to the 
purely interaction driven process (|^ for which the characteristic shape 





cosh J — cosh vf 
cosh J — cosh V 


, with,^ 



2.2005, n 


0.5413, J 


1 (16) 


can be analytically approximated na, as compared in Fig. The measured shape 


is similar to this approximation and as /3 increases, this characteristic shape seems to 
describe large 7 systems as well. 
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Figure 9. Comparison of the characteristic shapes for systems of various /3 and 7 at a 
condensate volume of about 10® masses. The shapes are formed by rescaling the width 
and height of all measured condensate sample shapes and only then averaging them 
(this avoids averaging artifacts one would get when interchanging these two steps). The 
fill colour inside the condensate shapes encodes the respective measured condensate 
width while the background colour around the shapes gives the critical density of the 
system. The shapes in the single-site condensate regime are plotted narrowed to give 
better distinction to extended shapes. 



Figure 10. Comparison of the rescaled condensate shape for /? = 1 ,7 = 0.2, M = 1000 
with the exact shape (16) derived in Ref. [T^ (blue line) for the interaction driven 
process Q. 


The phase boundary in our measurements is not a clear separation at /? = 1 between 
the two phases, as the predictions suggest, but it seems to be smeared out in favour of 
smooth condensates for small 7 and in favour of rectangular condensates for large 7 . 
In the respective regions the measured shapes combine characteristics of the individual 
phases. Once we have steep edges with smooth domes and on the other side we have 
slowly rising edges with a flat plateau on top. A comparison of the hnite-size shapes 
corresponding to different condensate sizes given in Fig. 11 (a-d) indicates that the 
transition line will sharpen with increasing condensate size, supporting that we are 
dealing with hnite-size effects. However, the prediction only tells us in the leading order 
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Figure 11. Finite-size effects of the condensate shape near the rectangular/smooth 
transition. The total amounts of masses in the system corresponding to the plotted 
rescaled shapes are with increasing width M = 10"^, M = 10®, M = 10® for 7 = 0.9 
and M = 10®, M = 10"^ M = 10® for 7 = 0.3 respectively. For /3 < 1, the smooth 
condensate edges disappear for increasing particle numbers M while they stay for 

/ 3 > 1 . 


which phase dominates for inhnitely large condensate volumes, but not whether there is 
a crossing where both phases are equally likely and at which volume this would occur. 
Therefore, the deviation of scaling exponents from the prediction could also imply that 
those systems have higher-order corrections to the transition. 

4 . 2 . Condensate width scaling 

Figure [T^ shows the estimated scaling exponents for the width-scaling behaviour of the 
smooth and rectangular case, with the measured values listed in Table [ 1 } One can see 
that for /d > 1 the measurements align well with the predictions, except for small 7 
values, where the deviations can be the result of hnite-size effects. There, zero-range 
interactions are weak and the condensates span many sites, which eventually becomes 
unfeasible to simulate for large condensate volumes. In the case jS < 1 and large 7 , 
the estimates are very close to the predicted values. Combining the scaling behaviour 
with the characteristic shapes, we see that for /? > 1.5 the condensate shape and scaling 
exponent is close to the behaviour of purely interaction-driven condensation, even for 7 
close to 1 . 

One region where the measurements deviate strongly from the predictions is for 
small 7 in the rectangular case. This region is also where the scaling behaviour is 
predicted to change dramatically at the phase boundary. It looks like the measurements 
of scaling exponents change in the same way as the predictions, but shifted to lower 
values. We have found no explanation so far why these systems behave this way. The 
obvious measurement outlier at position f3 = 0.9 and 7 = 0.1 can however be explained 
as an extension of the smooth scaling behaviour, which agrees with the observation that 
the condensate shape has no rectangular characteristics. 

The region at /d = 1 is of special interest, as there the condensate not only gradually 
changes from a smooth to a rectangular shape as 7 grows from 0 to 1 , but covers the 
whole scaling exponent range from 0 to 0.5 as well. In this region it also becomes 
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Figure 12. The measured scaling behaviour in the rectangular and smooth phases as 
a function of /?. In contrast to the phase diagram in Fig. the scaling exponent a is 
not shown in a ‘depth’ dimension, but as the y-value. Points for the same value of 7 are 
connected by broken lines. The continuous colored lines correspond to the predicted 
values for the scaling exponent a given by Eq. (14). The underlying condensate volumes 
are in the range of 3 x 10^ to 10^, depending on how extended the condensates and 
therefore how computationally expensive the simulations are. 


feasible to simulate the original dynamics and thus investigate dynamical properties of 
these systems as well. 

The peak condensate region fd < 1,7 > fd, with the transition line /9 = 7 is 
reproduced sufficiently well with problems at small values of 7 due to exceedingly large 
critical density as visible in the background colour code in Fig. Here we observe 
zero-range like condensates as the local interaction term is strongly dominant. 

5. Summary 

Based on the improved simulation method, especially in the rectangular phase, we were 
able to show for the pair-factorized steady state (PFSS) model ([^ that the three phases 
predicted in Refs. mm with peak-like, rectangular and smooth (bell-like) condensate 
shapes exist in about the regions they are expected in. By analyzing the characteristic 
shapes, we found that the hnite systems we considered show hnite-size effects that smear 
out the phase boundary between rectangular and smooth, bell-like shaped condensates. 

With the improved estimator for the condensate scaling behaviour we were able to 
conhrm that the scaling behaviour is clearly different from that of purely interaction- 
driven condensation processes [10] . Table demonstrates that the estimates for large 
regions of the parameter space are close to the expected scaling behavior. In the smooth 
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Table 1. Table of measured values of scaling exponents a (upper values) and their 
respective absolute deviations from the predicted values (lower values). For the referred 
theoretical predictions see Eq. (14). To estimate these scaling exponents we performed 
on the order of 10® Monte Carlo sweeps to take about 10® samples. 
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0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 

0.1 

0 . 00 ( 10 ) 

+.00 

0 . 00 ( 10 ) 

+ .00 

0 . 00 ( 10 ) 

+.00 

0 . 00 ( 10 ) 

+.00 

0 . 00 ( 10 ) 

+.00 

0 . 00 ( 10 ) 

+.00 

0 . 00 ( 10 ) 

+ .00 

0 . 00 ( 10 ) 

+ .00 

0 . 00 ( 10 ) 

+ .00 

0.2 

0 . 05586 ( 68 ) 

-.14677 

0 . 00 ( 10 ) 

+ .00 

0 . 00 ( 10 ) 

+ .00 

0 . 00 ( 10 ) 

+.00 

0 . 00 ( 10 ) 

+.00 

0 . 00 ( 10 ) 

+ .00 

0 . 00 ( 10 ) 

+ .00 

0 . 00 ( 10 ) 

+ .00 

0 . 00 ( 10 ) 

+ .00 

0.3 

0 . 01809 ( 98 ) 

-.14857 

0 . 0389 ( 33 ) 

-.0520 

0 . 00 ( 10 ) 

+ .00 

0 . 00 ( 10 ) 

+.00 

0 . 00 ( 10 ) 

+.00 

0 . 00 ( 10 ) 

+.00 

0 . 00 ( 10 ) 

+ .00 

0 . 00 ( 10 ) 

+ .00 

0 . 00 ( 10 ) 

+ .00 

0.4 

0 . 1551 ( 67 ) 

-.0757 

0 . 13103 ( 78 ) 

-.03564 

0 . 06497 ( 12 ) 

-.02594 

0 . 00 ( 10 ) 

+ .00 

0 . 00 ( 10 ) 

+.00 

0 . 00 ( 10 ) 

+ .00 

0 . 00 ( 10 ) 

+ .00 

0 . 00 ( 10 ) 

+ .00 

0 . 00 ( 10 ) 

+ .00 

0.5 

0 . 2260 ( 18 ) 

-.0597 

0 . 2017 ( 11 ) 

-.0291 

0 . 15376 ( 32 ) 

-.01291 

0 . 079277 ( 96 ) 

-.011632 

0 . 00 ( 10 ) 

+.00 

0 . 00 ( 10 ) 

+ .00 

0 . 00 ( 10 ) 

+ .00 

0 . 00 ( 10 ) 

+ .00 

0 . 00 ( 10 ) 

+ .00 

0.6 

0 . 2927 ( 11 ) 

-.0407 

0 . 26195 ( 75 ) 

-.02376 

0 . 21893 ( 29 ) 

-.01184 

0 . 16164 ( 12 ) 

-.00503 

0 . 087798 ( 21 ) 

-.003111 

0.00(10) 

+ .00 

0 . 00 ( 10 ) 

+ .00 

0 . 00 ( 10 ) 

+ .00 

0 . 00 ( 10 ) 

+ .00 

0.7 

0 . 34117 ( 87 ) 

-.03383 

0 . 31238 ( 27 ) 

-.02096 

0 . 27406 ( 27 ) 

-.01166 

0 . 22512 ( 11 ) 

-.00565 

0 . 165949 ( 43 ) 

-.000718 

0 . 089215 ( 94 ) 

-.001694 

0 . 00 ( 10 ) 

+ .00 

0 . 00 ( 10 ) 

+ .00 

0 . 00 ( 10 ) 

+ .00 

0.8 

0 . 38282 ( 67 ) 

-.02895 

0 . 35627 ( 38 ) 

-.01873 

0 . 32126 ( 31 ) 

-.01208 

0 . 27958 ( 12 ) 

-.00613 

0 . 227928 ( 52 ) 

-.002841 

0 . 168049 ( 28 ) 
+ .001382 

0 . 097067 ( 61 ) 

+.006158 

0 . 00 ( 10 ) 

+ .00 

0 . 00 ( 10 ) 

+ .00 

0.9 

0 . 4623 ( 10 ) 

+ .0178 

0 . 39033 ( 47 ) 

-.02143 

0 . 36312 ( 34 ) 

-.01188 

0 . 32649 ( 18 ) 

-.00684 

0 . 282072 ( 77 ) 

-.003642 

0 . 229113 ( 33 ) 

-.001656 

0 . 163916 ( 12 ) 

-.002750 

0 . 081862 ( 50 ) 

-.009047 

0 . 00 ( 10 ) 

+ .00 

13 1.0 

0 . 46611 ( 67 ) 

-.00758 

0 . 43925 ( 58 ) 

-.00520 

0 . 40296 ( 51 ) 

-.00880 

0 . 36059 ( 46 ) 

-.01441 

0 . 32191 ( 51 ) 

-.01142 

0 . 27902 ( 31 ) 

-.00670 

0 . 229285 ( 63 ) 

-.001484 

0 . 165404 ( 20 ) 

-.001262 

0 . 08919 ( 11 ) 

-.00172 

1.1 

0 . 47105 ( 59 ) 

-.00514 

0 . 44486 ( 49 ) 

-.00514 

0 . 41672 ( 43 ) 

-.00433 

0 . 38415 ( 37 ) 

-.00474 

0 . 34759 ( 32 ) 

-.00535 

0 . 30853 ( 27 ) 

-.00397 

0 . 26286 ( 22 ) 

-.00381 

0 . 21143 ( 18 ) 

-.00285 

0 . 15211 ( 17 ) 

-.00173 

1.2 

0 . 47340 ( 57 ) 

-.00487 

0 . 45060 ( 46 ) 

-.00394 

0 . 42709 ( 37 ) 

-.00148 

0 . 39673 ( 31 ) 

-.00327 

0 . 36528 ( 35 ) 

-.00314 

0 . 33210 ( 35 ) 

-.00123 

0 . 29166 ( 32 ) 

-.00245 

0 . 24856 ( 31 ) 

-.00144 

0 . 19675 ( 27 ) 

-.00325 

1.3 

0 . 47459 ( 55 ) 

-.00541 

0 . 45386 ( 44 ) 

-.00448 

0 . 43267 ( 36 ) 

-.00211 

0 . 40755 ( 27 ) 

-.00154 

0 . 37863 ( 24 ) 

-.00232 

0 . 34801 ( 20 ) 

-.00199 

0 . 31296 ( 22 ) 

-.00283 

0 . 27588 ( 20 ) 

-.00190 

0 . 23271 ( 18 ) 

-.00258 

1.4 

0 . 47652 ( 50 ) 

-.00497 

0 . 45820 ( 39 ) 

-.00334 

0 . 43626 ( 35 ) 

-.00374 

0 . 41410 ( 25 ) 

-.00257 

0 . 38885 ( 22 ) 

-.00245 

0 . 36133 ( 18 ) 

-.00231 

0 . 33114 ( 13 ) 

-.00219 

0 . 29712 ( 12 ) 

-.00288 

0 . 26078 ( 12 ) 

-.00237 

1.5 

0 . 47755 ( 48 ) 

-.00521 

0 . 45842 ( 39 ) 

-.00586 

0 . 44011 ( 33 ) 

-.00433 

0 . 41960 ( 25 ) 

-.00348 

0 . 39772 ( 19 ) 

-.00228 

0 . 37261 ( 17 ) 

-.00239 

0 . 34580 ( 12 ) 

-.00203 

0 . 316235 ( 100 ) 

-.001947 

0 . 28340 ( 11 ) 

-.00231 

1.6 

0 . 47679 ( 44 ) 

-.00708 

0 . 45910 ( 36 ) 

-.00757 

0 . 44299 ( 31 ) 

-.00529 

0 . 42405 ( 23 ) 

-.00453 

0 . 40435 ( 19 ) 

-.00306 

0 . 38183 ( 16 ) 

-.00278 

0 . 35761 ( 12 ) 

-.00239 

0 . 33100 ( 10 ) 

-.00233 

0 . 30196 ( 11 ) 

-.00239 

1.7 

0 . 47782 ( 46 ) 

-.00703 

0 . 46184 ( 35 ) 

-.00691 

0 . 44528 ( 29 ) 

-.00633 

0 . 42823 ( 22 ) 

-.00510 

0 . 40938 ( 19 ) 

-.00441 

0 . 38965 ( 13 ) 

-.00321 

0 . 36778 ( 12 ) 

-.00260 

0 . 344148 ( 95 ) 

-.002006 

0 . 31735 ( 11 ) 

-.00265 

1.8 

0 . 47948 ( 46 ) 

-.00624 

0 . 46234 ( 35 ) 

-.00824 

0 . 44750 ( 25 ) 

-.00705 

0 . 43121 ( 23 ) 

-.00629 

0 . 41431 ( 18 ) 

-.00504 

0 . 39633 ( 13 ) 

-.00367 

0 . 37680 ( 10 ) 

-.00251 

0 . 354548 ( 98 ) 

-.002594 

0 . 33038 ( 10 ) 

-.00295 

1.9 

0 . 47933 ( 45 ) 

-.00716 

0 . 46190 ( 34 ) 

-.01032 

0 . 44962 ( 27 ) 

-.00753 

0 . 43441 ( 22 ) 

-.00677 

0 . 41868 ( 18 ) 

-.00556 

0 . 40221 ( 12 ) 

-.00404 

0 . 38326 ( 12 ) 

-.00383 

0 . 363798 ( 92 ) 

-.002868 

0 . 341784 ( 99 ) 

-.003044 
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phase, the measured scaling exponents are within 2% of the analytically predicted values, 
which verifies that the assumptions and approximations in the analytical predictions of 
Ref. [I2] are justified. In fact, the prediction yields very good values of the scaling 
exponents even for finite systems provided no transition line is close. Systematic 
deviations of the transition line between the rectangular and smooth condensate shape 
regimes that we observed are largely due to finite-size effects. 
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